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We compare current-current correlators in lattice gauge calculations with correlators in different 
potential models, for a pseudoscalar charmonium in the quark-gluon plasma. An important ingre- 
dient in the evaluation of the current-current correlator in the potential model is the basic principle 
that out of the set of continuum states, only resonance states and Gamow states with lifetimes of 
C ' sufficient magnitudes can propagate as composite objects and can contribute to the current-current 

' correlator. When the contributions from the bound states and continuum states are properly treated, 

the potential model current-current correlators obtained with the potential proposed in Ref. [11] 
are consistent with the lattice gauge correlators. The proposed potential model thus gains support 
<*"] ' to be a useful tool to complement lattice gauge calculations for the study of QQ states at high 

, temperatures. 



PACS numbers: 25.75.-q 25.75.Dw 



(N ' I. INTRODUCTION 

> 

Potential models have often been used to describe bound states of quark-antiquark pairs. The basic idea is that a 
^-j- . short range attractive color-Coulomb interaction with a long-range confining interaction provides an adequate account 
of the interaction between a quark and antiquark. While the non-relativistic potential model was originally introduced 
t-H , for heavy quarkonium systems, relativistic and non-relativistic quark models using constituent quark masses have been 
^»<0 ' used to describe mesons with one or two light quarks as quark-antiquark bound states [H, 0, 0, 0, 0] ■ 

Potential models have also been used to study heavy quarkonium bound states at high temperatures [H, 0, 0, 
0, EH, El, EH, EH, EH, EH, EH, EH, EH, EH]- At temperatures above the phase transition temperature, the potential 
Oh! between a quark and an antiquark is subject to screening [f|. Heavy quarkonium states become unbound in the 
q. potential as temperature rises. Potential models provide beneficial complements to lattice gauge calculations as 
' potential models allow simple and intuitive evaluations of many quantities, some of which may be beyond the scope 
C"| , of present-day lattice gauge calculations. The central question in the potential model has been focused on the 
characteristics and the temperature dependencies of the screening potential as determined by lattice gauge calculations 

0, 0, 0, EH, EH, El, EH, EI EH, EH, E3, El, El . 

■ There is recently a serious theoretical question whether it is appropriate to apply a potential model to study 



heavy quarkonia at high temperatures |2Ct [2l|. On the one hand, lattice gauge spectral function analyses have 
been carried out to investigate the stability of heavy quarkonia at high temperatures [H, [5H, HH, HE HH| ■ On the 
other hand, independent lattice gauge calculations [27], [28( have been performed to evaluate various thermodynamical 
quantities, such as the free energy -Fi(r, T) and the internal energy C/i(r, T), for a color-singlet static Q-Q pair at 
various separations r and temperatures T. ft remains an important open theoretical question how a quark-antiquark 
potential can be extracted from these thermodynamical quantities. Various potential models have been proposed, 
utilizing Fi @, 0, El, U± 9l. fiol. 16l . EU, or a linear combination of Fi and U\ with coefficients that depend on the 
equation of state [ll|, El, EH, Ei|- Although the latter linear-combination model has been justified by theoretical 
arguments and leads to dissociation temperatures that are consistent with lattice gauge spectral function dissociation 
temperatures [ll|, El, EH, EH] , it remains necessary to confront the model with other lattice gauge data to assess the 
degree of its usefulness. 

The spectral function is related to the meson current-current correlator by a generalized Laplace transform. In prin- 
ciple, they carry equivalent information on the composite system. One would expect intuitively that the consistency of 
the lattice gauge spectral function dissociation temperatures with the potential model dissociation temperatures using 



the potential of would lead to consistency of the lattice gauge correlators with the corresponding potential model 



correlators. Recent works in Refs. [20j, |21| however make the contrary claim that the meson correlators obtained with 
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many different types of potential models are not consistent with lattice gauge correlators, and consequently potential 
models cannot describe heavy quarkonia above T c . 

The failure of the potential model correlators to reproduce lattice gauge correlators in Rcfs. [2(| HU for all cases 
suggests that the lack of agreement may not be due to the potential models themselves but to the method of evaluating 
the meson correlators in the potential model. In the work of Ref. [2(| HH , continuum states arising from a free fermion 
Q and Q pair in a fermion gas contribute to the correlator. However, to be able to propagate as a composite meson 
so as to contribute to the correlator, the quark and the antiquark must be temporally and spatially correlated to be 
a composite object with a sufficiently long lifetime. Continuum states in the free quark and antiquark gas may not 
have sufficient temporal and spatial correlations to be composite objects for such a propa gati on. 

While we raise questions on the method of evaluating the potential model correlators in |2dl . [2l| , we wish to present 
what we view as a proper treatment of the current-current correlator in the potential model. We wish to point out 
the basic principle that out of the continuum states, only resonance states and Gamow states [29j with lifetimes of 
sufficient magnitude can propagate as composite objects and contribute to the meson current-current correlator. With 
the simple example of the pseudoscalar correlator, we shall show in this paper that when both the bound state and the 
continuum state contributions are properly treated, the potential model of [Til ] using a linear combination of i*i(r, T) 
and U\{r,T) yields correlators consistent with lattice gauge correlators, while the F\ potential and the U\ potential 
lead to deviations. Our results indicate consistency of the potential model of with both the lattice gauge spectral 
function analyses and the lattice gauge correlator analyses. The potential model of [Til ] thus gains support to be a 
useful tool to complement lattice gauge calculations for the study of QQ states at high temperatures. 

In Section II, we review the basic assumptions in the treatment of the continuum states in Refs. [2(| We 
review in Section III the relationship between the current-current correlator and the quarkonium wave functions. In 
Section IV, we show how the resonance states and Gamow states are characterized in the potential model. In Section 
V, the potential model correlator is expressed as a sum over contributions from meson bound and continuum states 
with lifetimes greater than a certain limit. In Section VI, we discuss the relation between the Q-Q potential and 
lattice gauge thermodynamic quantities. In Section VII, we evaluate the potential model correlators and find that the 
correlators obtained with the potential of [ll[ have features similar to those of lattice gauge correlators. We have thus 
demonstrated the consistency of the potential model correlators with lattice gauge correlators. In Section VIII, we 
discuss the implications of the potential model analysis on the assume default spectrum in the lattice gauge spectral 
function analysis. In Section IX we present our conclusions. 



II. MESON CORRELATORS AND CONTINUUM STATES 



The meson (current-current) correlators G(r, X) is a function of the Euclidean time r and the spatial coordinate 
X defined by 

G(t,X) = (J m (t,X)jI(0,O)), (1) 

where Jm(t, X) = g(r, X)Fm<;(t, X) and Tm = 1, 7 5 , 7 M , 7 5 7' 1 are the operators appropriate for scalar, pseudoscalar, 
vector, and axial-vector mesons respectively. It is the probability amplitude for creating a meson M at space-time 
point (0,0), propagating the meson to (r, X), and destroying the meson at (r, X). Specializing to the case of the 
meson momentum P equal to 0, the meson correlator depends then on the meson mass spectrum er(a;, T). 

In their test of the potential model, the authors of 0, [2l[ assume that the meson mass spectrum in the potential 
model is given by 

a(<j,T) = ]T 2M. l E l S(u; 2 - M?) + -^6{u - s )(a H + (2) 

i 

where Mi(T) is a bound state meson mass calculated in the potential model, Fi(T) is the corresponding magnitude of 
the wave function at the origin. The second term with the step function represents continuum meson states obtained 
by assuming that the Q and Q are free fermions in a fermion gas [3C| . The quantity sq{T) is the continuum threshold, 
au and bn are constants that depend on the meson type. The potential model meson correlator G(t) is then obtained 
by folding the meson mass spectrum <t{uj,T) with the propagating kernel K{t,uj,T), 

G(t) = J o!ujo-(lu,T)K(t,lu,T), (3) 

where K(t, ui, T) is given by 



(4) 
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The lattice gauge correlator G(t) is represented relative to the "reconstructed" correlator G recon (r) calculated with 
the meson mass spectrum at T — 0, [3l[ . The lattice gauge pseudoscalar correlators obtained in [24[ are shown 
in Fig. la. The potential model correlators G(t)/G YCCOVl (t), obtained in 0, HH with a screening potential (Fig. lb) 
or with the potential of Ref. [HI, [H (Fig. lc), are found to be qualitatively very different from the lattice gauge 
correlators of Fig. la. Significant differences occur for both charmonia and bottomia, at all temperatures above T c , 
and for many fundamentally different potentials. The authors of Ref. 0, [2l[ then draw the conclusion that the 
potential model is not a good description for heavy quarkonia above T c . 




FIG. 1: Fig. la shows the pseudoscalar charmonium correlators G{t) / G TCCO u{t) obtained in lattice gauge calculations [24l.l3ll[. 
Fig. lb shows the potential model correlators for the pseudoscalar charmonium obtained in Ref. [2(| with the screened Cornell 
potential. Fig. lc gives the potential model correlators for the pseudoscalar charmonium obtained in Ref. [2lJ with the potential 
of Refs. |H Cl- 
in the evaluations of the potential model correlators in [20l . l2l| , contributions from continuum states represented by 
the second term in Eq. ([2]) have been included. While the spectrum in the continuum states is continuous as in Eq. 
([2]), not all continuum states represented by the step function QQ free-fermion gas states in the second term in Eq. 
@ possess the proper characteristics to propagate as a composite meson so as to contribute to the meson correlator 
G(r). To be able to propagate as a composite meson from time to r so as to contribute to the correlator G(r), 
the quark and antiquark must be temporally and spatially correlated. The composite object must have a composite 
lifetime sufficiently long compared to the propagation time r. Continuum states in the free quark and antiquark gas 
may not have sufficient temporal and spatial correlations between Q and Q to qualify as being composite objects for 
such a propagation. 

A similar question poses itself in low-energy nuclear physics in the evaluation of the density of continuum states of 
a composite object formed by a particle (neutron, say) and a nucleus (represented by a potential well). States for the 
free particle in a (neutron) fermion gas do not have relative temporal and spatial correlations between the particle 
and the nucleus to qualify as being a composite object formed by the particle and the nucleus in the continuum. It is 
for this well-known reason that to get the level density of continuum states of a composite object formed by a particle 
and a nucleus, it is necessary to subtract the level density of free gas states from the total density of continuum states 
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[32l . [33l . l34l . [35l . l36l . [3^ . In the analogous problem of level density of continuum states of a composite QQ meson, one 
can carry out the free fermion gas states subtraction or alternatively use only resonance states and Gamow states 
as contributors to the continuum level density [H, [H, [H, HI, H3- Therefore, for our case of meson correlators 
for a composite QQ meson, the basic principle is that the proper continuum meson states that can contribute to the 
meson correlator G(t) should be limited to meson resonance states and Gamow states which have composite object 
lifetime long compared to r. By focusing our attention on the pseudoscalar charmonium as an example, we would 
like to demonstrate that the potential model with the correct set of bound and continuum Gamow states can describe 
lattice gauge correlators. 



III. THE MESON CORRELATOR AND MESON INTERNAL WAVE FUNCTIONS 

We start with the current-current correlator G(t, X) = ( Ja/(t, X) J* M (0, 0)) of Eq. |T]) and restrict ourselves to the 
consideration of the pseudoscalar charmonium. The current operator is just Jm{t, X) = q(r, X)75g(r, X). The field 
operator has both annihilation and creation parts, g(r, X) = q^ + \r, X) + q^' (r, X). Hence, the current operator 
can be decomposed into 

Jm(t, X) - (r, X) l5 q i+) (t, X) + g« (r, X) l5 q^ (r, X) + g<"> (r, X) 75 g« (r, X) + g« (r, X) W~> (r, X )- (5) 

We consider the propagation of the pseudoscalar meson M from the initial (Euclidean) space-time coordinate (0, 0) 
to the final space-time coordinate (t, X) . It is often convenient to study 

G(r) = y"dXG(r,X)= J dX( J m (t, X) 4,(0, 0)). (6) 

The temporal behavior of the correlator G(r) provides information on the meson spectrum. 

For the extraction of information on the correlator we focus only on the particle-antiparticle creation operator 
portion q( + *> (t, X)75g( _ )(T, X) of Eq. The 75 matrix singles out the dominant lower component of the operator 
<P'{t, X) for the antiquarks. We rename that dominant component by using the notation 



Jm(t,X) = 1pl(xi)lj)2(X2] 



(7) 

X\ — X2 — (T,X) 



where particle 1 represents Q and particle 2 represents Q, ipi(xi) and ^2(^2) are the dominant Q and Q fields. The 
meson current-current correlator of Eq. |T]) becomes 



G(r,X) = ( Mxi)M^)4^'i)4(^) ) 



(8) 

x\— X2 — (r;X);a:^— x' 2 — (0,O) 



The above current Jm{t, X) is a local current where (r, X) = X is the CM coordinate X of the Q-Q pair with 
X = {x\ + X2)/2, and the relative coordinate, x = (x\ — X2), is restricted to be zero. A general current containing a 
more general relative coordinate x = (xx — X2) is 

Jm(X,x) = MX + x/2)MX - x/2), (9) 

with the more general Green's function 

G(x l ,x 2 ;x[,x' 2 ) = {^ 1 {x l )^ 2 {x 2 )^\{x' 1 )^ 2 {x l 2 )). (10) 

We would like to make a transformation of the coordinate system from the separate particle coordinates {x\,x 2 ) to 
(X, x) and introduce the field operator ^m(X) to represent the field of a composite meson which has the internal 
relative motion between the quark and the antiquark described by an internal wave function ipM- 
Following [38[, we define the field operator ^m(X) for CM motion by 

9 M (X)\X) = * M (X)\X), (11) 
with the wave function ^m(X) satisfying Eq. (2.23) of [38| 

{P 2 -M 2 }^ M (X) = 0, (12) 
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where P — pi + P2 is the total momentum of Q and Q, pi is the momentum of of Q, and p 2 is the momentum of Q. 
For each composite state of mass M, there is a wave function iPm(,x±) for relative motion satisfying Eq. (2.24) of [38[ 

{p 2 - ^(x ± ,M)+b 2 (M 2 ,m 2 ,m 2 )}^ M (x ± ) = 0, (13) 

where p = (pi — p 2 )/2 is the relative momentum of Q and Q, $(xj_,M) is the two-body interaction potential that 
acts between Q and Q, and 

PP ' X HAS 

xx = x -p^-, (14) 

P-x± = 0. (15) 

The composite particle mass M is the constant of the separation of variables and is related to the eigenvalue b 2 of the 
Schrodinger-type equation for relative motion, Eq. (|13[) . by 



6 2 (M 2 , m 2 , m 2 ) = ^{M A - 2M 2 (m 2 + m 2 ) + (m 2 - m 2 ) 2 }. (16) 

We decompose the current field operator as a sum of composite meson operators ^ n (X) with coefficients ij)n(x±), 

Jm(X, x±) = Mx±)^n(X), (17) 

n 

in which the summation n includes a sum over bound and continuum states. As the field operators of composite 
objects of different types and energies produce orthogonal states, we get 

G(x 1 ,x 2 ;x[,x' 2 ) = G(X,X';x ±1 x' ± ) =^ n (x ± )i>*(x' 1 _)(^ n (X)4>l(X')}. (18) 

n 

The Green's function for a composite meson of mass M n is 

g n (X,X') = {* n (X)& n (X% (19) 

we have therefore 

G(X,X';x x ,x' 1 _) = J2M^Wn(x'x)9n(X,X'). (20) 

The advantage of the separation of the correlation function in terms of the relative degrees of freedom {x± , x' ± } and 
the composite particle CM degrees of freedom {X, X'} is that the Green's function g n (X,X') corresponds to that of 
a free single particle (meson) of mass M n in a thermal bath of the plasma We consider X 1 = and X = (r, X) 
with < t < (3 — l/T and take the spatial Fourier transform of <7„(rX, 0) with respect to X, 

g n (r,P)= y"dX e - lPX 5ti (rX,0). (21) 

The correlator of a free particle with < r < (3 is given by [see e.g. the equation before (3.4) of Ref. [39j |] 

. . _ f dPon PanT A(P 7 P 0n ) 

5 " (T ' P) -i 2 m e i_exp{-/3(P 0n -M)}' [ ' 

where [i is the chemical potential for the composite particle in the medium and A(P,P ) is the spectral function for 
a free composite particle given by [see e.g. equation after (3-8b) of [3!| for the non-relativistic case], 



A(P, P 0n ) = 2TrS{P 0n - v/P 2 + M2). (23) 
The Green's function for the composite meson in the state n with P = is then 

~g n {r) = 3n(r,P = 0) = fdX g n (rX) = - ^f'fjf n (24) 



G 



The corresponding correlator G(r) in Eq. (O and (f2T))) is then 



G(r; x±,x' ± ) = \ £ ip n (xj_ ,M n )ip* (x' ± , M n ) j ^ (25) 

In the center-of-mass system of the composite particle at rest, x± = (0,r), and we have 

G(r; r.-O = J E Mr, M„)«r', M„) , _ (26) 

In the above equation, the propagator kernel for a composite particle in state n with a mass M n can be re-written as 

- n = exp(-M ra r) = cxp{-M n (r-/3/2)+/3/z/2} 

5nlTj 1 - exp{-/3(M„ — /it)} 2sinh{/3(M n -/i)/2} ' 1 j 

In lattice gauge calculations, one chooses to work with the case of [i = and imposes the periodic boundary condi- 
tion, G(r) \ T= o = G(t) \t=i3- The periodic boundary condition can be satisfied by including not only the exponentially 
decreasing solution e~ M " r but also the exponentially increasing solution e M " r in Eq. (j22|) so that the propagating 
kernel becomes, 

K( M T s _ exp{-M»(r - (3/2)} + cxp{M n (r - 0/2)} _ cosh{M ra (r - (3/2)} 
l ' >~ 2sinh{/3Af„/2} ~ sinh{/3M„/2} ' 1 > 

which satisfies the periodic boundary condition, K(0, M n ,T) = K((3, M n ,T). We shall consider this case of ji = 
with the periodic boundary condition and use the above propagating kernel, in order to compare the potential model 
correlators with the lattice gauge correlators. We shall however compare results only for r less than and up to (3/2 
before the onset of the dominance of the exponentially increasing evolution. For the case considered in lattice gauge 
calculations, where the relative coordinate of the Q-Q pair is set to zero, we have r = r' = and thus 



G(r)=G(r;r,r') 



- V K(t, M n , T)^„(r,M n )C(r', M n ) 

7 ^ ^ 



(29) 



IV. GAMOW STATES AND RESONANCE STATES IN THE CONTINUUM 



We would like to use Eq. (|29|) to evaluate correlators for the pseudoscalar charmonium in the potential model, for 
comparison with pseudoscalar correlators obtained in lattice gauge calculations. As the mass of a charm quark is 
quite large, we shall restrict ourselves to non-relativistic considerations. For simplicity, we shall also ignore spin. The 
summation over n can be separated into a sum over bound states b and a sum over continuum states k. The proper 
treatment of the contributions from the continuum states is important in the evaluation of the correlator. 

As we remarked earlier, the meson correlator G(r, X) is the probability amplitude for creating a meson M at 
space-time point (0, 0), propagating the meson to (r, X), and destroying the meson at (t, X). To be able to propagate 
as a composite object from time to r, the quark and antiquark pair must be spatially and temporally correlated. 
The composite object must have a lifetime long compared to the propagation time r. Using considerations similar 
to those in the density of states of a composite object in nuclear physics [H, [H, H3, HH, [H, |37|, we are well advised 
that the proper continuum states that can contribute to the meson correlator G(r, X) should be resonance states and 
Gamow states with lifetimes long compared to r. 

We need information on the lifetime of a continuum state k in the potential model. For that purpose, we can 
calculate the phase shift <5(ek) as a function of the continuum state energy ek = k 2 /2/^ m where \i m is the reduced 
mass, mo/2. Kno wing the phase shift as a function of energy, we can determine the delay time, as pointed out by 
Wigner [iol. [4ll. 1421. |43|| . To obtain such a relationship, we consider the wave packet with momenta centered around k 
with energy ek to travel from r = to R with R along k, as a function of t. The peak of the wave packet arises from 
the interference of two waves with energy ek and ek' • They must interfere constructively. Constructive interference is 
possible when the phase difference of the two wave functions at ek and ek' is zero. This condition for the constructive 
interference can be written explicitly as 

k R k' R + 6(e k ) - <5(e k <) - (e k t - = 0. (30) 

Solving for t and taking the limit k' approaches k, we obtain 

de^/dk 9ek ' ^ ^ 
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which indicates that the passage of the wave packet with momentum centered at k from the origin to R is delayed 
by a delay time [io| 



(delay time) = 



d6(e k ) 
de k 



(32) 



A negative delay time, with a negative cW(ek) / de k , represents the flying apart of the Q and Q in advance of their coa- 
lescence approach and by causality cannot represent a composite object. Continuum states with negative cW(ek)/<9ek 
should not be included as contributing to the meson correlator for the propagation of the composite quarkonium. 

A system with a positive delay time can be interpreted as possessing a finite lifetime. What fraction of the delay time 
should be assigned to the lifetime of this composite object residing as a wave packet centered around the continuum 
state k? To answer such a question, we seek the help of the case of a sharp resonance for which the answer can be 
readily obtained. 

If the state at k is a sharp resonance, the phase shift <5k(e) in the neighborhood of the resonance k with a small 
width Tk is given by 



tan <5k (e) = 



r k /2 



e k - e 

Then the derivative of the phase shift with respect to the continuum energy is 

94(e) T k /2 



de 



(e k - e) 2 + [T k /2] 



2 ' 



At the resonance e = e k , we have 



dS k (e) 



dc 





T k " 




2 



(delay time). 



(33) 



(34) 



(35) 



Thus, by examining the case of a sharp resonance state, we find that half of the delay time can be assigned to 
the lifetime, l/r k , of the continuum state. We shall therefore assume that the assignment of the fraction of one-half 
of the delay time as the composite particle lifetime is a reasonable concept, for all states with positive <9<5 k (e k )/<9e k . 
Thus, we assign the width Tk associated with a state k with a time delay as 



r k = 



<9<S k (ek)/de k ' 



(36) 



A well-defined potential resonance in the continuum with momentum k and energy e k occurs when <9(5 k (e k )/<9e k > 
and <5 k (e k ) = (2n+l)7r/2, where n is an integer. States at energies with <5 k (e k ) = (2n+l)7r/2 but with <9<5 k (e k ) / 'de k < 
are echos and not resonances j44[. For these reasons, if the S- wave potential does not possess a barrier to trap the 
waves in the interior, there will be no S-wave resonances j4l| . Nevertheless at plasma temperatures above ~1.6T C (as 
we shall see in subsequent sections), there are S- wave continuum states with positive time delays and lifetimes and 
can be represented by Gamow states with various widths. They will be used in our subsequent calculations of the 
potential model correlators. 



V. MESON CORRELATOR IN THE POTENTIAL MODEL 



Upon limiting our attention in this manuscript to the pseudoscalar charmonium with L=0, we note that the quark- 
antiquark potential itself does not possess a potential barrier, and thus there arc no sharp 5-wave resonances. There 
are however continuum states with positive time delays at T >~ 1.6T C , as we shall see in Section VII. These continuum 
states can be represented by Gamow states with finite lifetimes and widths. 

With composite particle states consisting of bound states and Gamow states, the summation over n in Eq. (|29|) 
for the evaluation of the correlator in the potential model can be separated into a sum over bound states b and an 
integral over continuum Gamow states k, 

M*, M n )r n (r', M n )f(M n ) / dwj^ S(u - M 6 )^(r, M b ))^(r', M b )f{uj) 

n J b 

+ dej dk^ k (e,r k )0(5 ( 5k/aek)^k(r)Vk(r')/(e) (37) 
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where the width of the Gamow state is 2/[<9<5k(ek)/<9e k )] and the function -D k (e,r k ) represents the spreading of 
the distribution of the Gamow states k in continuum energy e due to the presence of its delay time and a width r k . 
One can choose different forms of the distribution function -D k (e,Fk) that has a peak at £k with a width Tk, such 
as a Gaussian or a Breit-Wigner distribution in e. The results in the correlators would not be sensitive to the form 
of the distribution function -Dk(e,Tk). For convenience, we shall choose to represent -Dk(e,Tk) by the Breit-Wigner 
distribution 

. . CLfc Tk/2 . . 

Dk{e,T k ) — — , (38) 

7r (e - e k r + ( r k/2) 2 



where 



2-7T 

a , = (39) 

7r + 2arctan(2e k /r k )' k ; 



and J °° deD k (e, r k ) = 1. Eq. |37|) becomes 

X] M r > M nWn(r\ M n )f(M n ) -» fduJ2 S (u- Mb)M*, M 6 ))V$(r', M 5 )/(c 







„ a/c (r k /2)6»(a<5 k /ae k ) , 
de / rfk V ( £ - £k ) 2 + (r k /2)^ k(r) ^ (r )/(e) (40) 



We normalize the wave function of the bound state Mb according to 

dr|^(r J M 6 )| 2 = l. (41) 
The wave function in the continuum ?/» k (r) is normalized according to 

dr ^k(r) V k '(r)=5(k-k'), (42) 



and it behaves asymptotically as 



, , . expjik • r} , , . 



It is not sufficient to limit the width r k to be positive. To be able to propagate temporally as a composite particle 
from to r, the composite object must have a lifetime exceeding a certain limit h/T max . It is therefore necessary to 
limit the contributions in the integral in Eq. (|40| further by an additional step function ^(rmax — r k ). The correlator 
is therefore 



1 cosh{M 6 (r- 3/2)} , s , 



{$mh{M e 8/2} 

de I d^ ^ r k ) ca*{M{T-W)} 

J it (e - e k ) 2 + [r k /2] 2 sinh{M e( 3/2} W j 1 ; 

We need to specify the maximum width limit r max . The width limit T max depends on the time scale in the 
propagation of the composite meson. In lattice gauge calculations, the composite object propagation time varies from 
r = to t = 8/2 with a maximum of r = 3/2. A propagation time r greater than 3/2 will lead to the unphysical 
region where the meson probability amplitude grows predominantly exponentially with time (see Eq. (|28|) ). With this 
maximum r = 3/2 in the correlator measurement, the composite object lifetime riif e needs to be greater than 3/2 
and the minimum life time r m ; n nf c is given by 

r min nfc(T) = 3/2 = 1/2T. (45) 

This arises because a composite object with a lifetime Tuf e less than r m i n iif e (r) will not be able to remain a composite 
object as it propagates from r = to r = 8/2. For the continuum state at k, the above requirement, that the 
composite object lifetime riif e = 1/T k must be greater than or equal to r min nf c (T) = 1/2T, leads to the condition for 
the maximum width r max = l/r m i n uf e of the continuum states detectable by the correlator measurement, 

2T. (46) 



nun 



lifc(T) 8/2 
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For the case considered for the lattice gauge calculations where the relative coordinate x = {x\ — x 2 ) of the Q-Q 
pair is set to zero (Eq. ©), we have r = r' = and thus we need to evaluate G(t; rr') at r = r' = 0. We define the 
quantity K, as the the ratio of the absolute square of the amplitude of the continuum wave function at the origin to 
the absolute square of its amplitude at infinity [45[ 



/C(k) 



l'0k(O)| 5 



|V>k(r ->■ oo)| 2 ' 



(47) 



which can be evaluated numerically for the potential in question. Using Eq. (|43[) . the continuum wave function at the 
origin is therefore approximately 



|^ k (0)| 2 = /C(k)/(27T) 3 



(48) 



The correlator becomes 

G(r) = G(r; rr'; 

+ j / 



cosh{M h (r - /3/2)} 
sinh{M b /3/2} 



r=r'=0 " b 

a k (T k /2)6(dS k /de k )e(T max - T k ) cosh{M e (r - 0/2)} £(k) 



dk- 



(e - e k ) 2 + [r k /2] 5 



sinh{M £ /3/2} (2tt) 3 



(49) 



VI. RELATION BETWEEN THE Q-Q POTENTIAL AND LATTICE GAUGE THERMODYNAMIC 

QUANTITIES 

In the potential model, the most important physical quantity is the Q-Q potential between the quark Q and the 
antiquark Q in a color-singlet state. Previous works in the potential model use the color-singlet free energy Fi @,[II,E3 
or the color-singlet internal energy Ui 0, E3, EB, EH obtained in lattice gauge calculations as the color-singlet Q-Q 
potential without rigorous theoretical justifications. The internal energy U\ is significantly deeper and spatially more 
extended than the free energy F\. Conclusions will be quite different if one uses the free energy F\ or the linear 
combination of F\ and U\ defined by Eq. |58|) below as the Q-Q potential. 

If one constructs the Schrodinger equation for the color-singlet Q and Q, the Q-Q potential U^),(r,T) in the 

_ _ WW 

Schrodinger equation contains those interactions that act on Q and Q, when the Q and Q are separated by r at 
temperature T and the medium particles have re-arranged themselves self-consistently. As shown theoretically in 
detail in for lattice gauge theory, this potential is given by 

U$± (r , T) = Ux (r , T) - [U g (r, T) - U g0 ] , (50) 

where Ui(r,T) is the color-singlet internal energy, U g (r,T) and U g o are gluon internal energy in the presence and 
absence of the color-singlet Q and Q pair, respectively. This is a rather general result when screening occurs, as a 
similar relationship exists between the total internal energy and the Q-Q potential in the analogous case of Dcbyc 
screening [l3j ]. 

We proposed earlier a method to determine the gluon energy in Eq. (|50[) in terms of the gluon entropy by making use 
of the equation of state of the quark-gluon plasma obtained in an independent lattice gauge calculation in quenched 
QCD [ll], Qjl, EH- The equation of state of the medium provides a relationship between the QGP internal energy 
density e g and the QGP entropy density a, 

e g =Ta - p. (51) 

Thus, by expressing p as (3p/e s )(e g /3) with the ratio a(T) = 3p/e g given by the known equation of state in quenched 
QCD, the plasma internal energy density e g in quenched QCD is related to the entropy density To by 

To. (52) 



3 + a(T) 



This is just 



dV 3 + a(T) dV 



J dV T(a-a a + a ), (53) 
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where a® is the entropy density in the absence of Q and Q. Noting that the entropy of the medium for the color-singlet 
Q-Q pair is TSi = J dV T(a — a ) and U g o is related to J dV Ter , the above equation leads to 

dV 3 + a{T) dV [ ' 

The plasma internal energy difference integrated over the volume is therefore given by 

UW(T,T) - U g0 (T) = 3+ l (T f S 1 (r,T). (55) 

But TSi(r,T) has already been obtained as TSi(r,T) — U\(v,T) — F\(r,T) from the lattice gauge calculations [27| . 
The plasma internal energy difference, U g (r, T) — U g o, is therefore equal to 

C/ W (r , T) - U g0 = 3+ [ {T) Pi (r, T) - F 1 (r, T)], (56) 

The Q-Q potential, l/j^lfr, T), as determined from Eq. (TSDl) by subtracting the above plasma internal energy difference 
from U\{y,T), is then a linear combination of F\{v,T) and U\{y,T) given by [11| 

U { Q 1] Q (r,T) = J^rf^T) + j^U^T). (57) 



We can rewrite the above as 

W/. (y T\\ = T, 

QQ 



)) ee U%l(r,T) = f F (T)F 1 (r 1 T) + f u (T)U 1 (r,T) 1 (58) 



where for brevity of notation we have renamed Ui,L(r,T) as Wi(r,T) and we have defined the coefficients 

WW 



and 



f F (T) = (59) 
3 + a(T) 

fu(T) = ^^r. (60) 
3 + a(T) 

To determine a(T), we use the equation of state of Boyd et al. [46[ for quenched QCD. The values of a(T), f F (T), 
and fu(T) are given in Fig. (1) of Ref. [ll[. The W\ (r, T) potential is approximately Fi(r 7 T) near T c and is 
3Fi(r,T)/4 + [/i(r,T)_/4 for high temperatures [HI, [13, [3 . 

To evaluate the Q-Q potential Wi(r,T), we use the free energy Fi(r, T) and the internal energy Ui(r, T) obtained 
by Kaczmarek et al. in quenched QCD [27| for which -Fi(r, T) and U\(r, T) can be parametrized in terms of a screened 
Coulomb potential, 

Fx r,T) = C F (T) ^Jl , (61) 

3 r 

and 

U 1 (v,T) = C u (T)--^- r . (62) 

The parameters Cf,u{T), a F .u(T), and /i Ft u(T) arc shown in Figs. 2 and 3 of Ref. [TTI | . 
The non-relativistic Schrodinger equation for the QQ states in the Wi(r,T) potential is 

+ W 1 {v 1 T)+m Q +m Q X i/i(t,T) = M(T)i/>(r,T). (63) 



2/i m 

This equation can be re-written as 



+ Wi(t,T)- VFi(r^oo,T) \^>(r,T) = {M(T) - m Q - m Q - Wi{r -»■ oo, T)}^(r, T) 



2/' 



e(T)^(r,T). (64) 
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With the V^i(r, T) potential as given by Eq. (|58p. (|61[) . and (JHU), the mass of the composite meson and the eigencncrgy 
e(T) of the above equation are related by 

M{T)=m Q +m Q + e(T) + C(T), (65) 
where C(T) is the asymptotic value of the Wx(r — > oo,T)[= U^(r — > oo,T)] obtained in lattice gauge calculations, 

C(T) = Wi(r - oo, T) = f F (T)C F (T) + fu(T)Cu(T). (66) 

For simplicity, we ignore spin and we study the S-wave charmonium state which splits into J/ip and r\ c when the 
spin-spin interaction is included. Using the Wi(r, T) potential and m c = 1.41GeV, we calculate the bound L=0 
charmonium state energy e(T). The mass M(T) of the bound S state can then be obtained from e(T) using Eq. (|6"5"|) 
and C'(T) = Wi(r -» oo,T). We list C(T), e(T), and M(T) for the L=0 charmonium states in Table I. We also list 
M(T = 0) = 3.064 GcV given by M(S state, T = 0) = 3M{J/ip,T = 0)/4 + M(r? c ,T = 0)/4 where M(J/ip,T = 0) 
and M(rj c , T — 0) are experimental masses. 



Table I. The quantities C(T), e(T), M(T) (all in GeV) for the bound L=0 charmonium state, 
calculated with the Wi(r,T) potential of Eq. in quenched QCD. 



T 


0. 


1.13T C 


1.18T C 


1.25T C 


1.4T C 


1.6T C 


1.95T C 


2.60T C 


C(T) 




0.2962 


0.2710 


0.2476 


0.2218 


0.1991 


0.1620 


0.1094 


Bound state energy e(T) 




-0.0340 


-0.02078 


-0.0105 


-0.0036 


-0.00019 


unbound 


unbound 


Mass M(T) 


3.064 


3.082 


3.070 


3.057 


3.038 


3.019 



As one observes, the bound S state masses obtained in the potential model with the potential of Eq. (|58|) are nearly 
the same as the S-wave charmonium mass at T = 0. They change only very slightly with temperature until the 
S-wave charmonium dissociates at ~1.6T C . The dominant errors in the mass value M(T) comes from the statistical 
errors in the asymptotic values Cf,u(T) of F\(r, T) and Ui(r, T) in [27j . They leads to uncertainties of about ±0.7% 
in the values of M(T). 

Table II. Spontaneous dissociation temperatures calculated with the W\{r,T) potential of Eq. (|58[) . 
the F\{r, T) potential, and the U\{t,T) potential, in quenched QCD. 





Dissociation temperatures in Quenched QCD 


States 


Spectral Analysis 


Wi Potential 


Fi Potential 


U\ Potential 


J/ip,Vc 


1.62-1. 7072 


1.62 T c 


1.40 T c 


2.60 T c 


Xc 


below LIT,? 


unbound in QGP 


unbound in QGP 


1.18 T c 


V,Tfc 




unbound in QGP 


unbound in QGP 


1.23 T c 






4.1 T c 


3.5 T c 


-5.0T C 


Xb 


1.15-1.54TJ 


1.18 T c 


1.10 T c 


1.73 T c 






1.38 T c 


1.19 T c 


2.28 T c 



fRcf.[22], ^Ref.[24j, »Ref.[25, 26] 



In lattice gauge spectral function analyses, the positions of the bound S states do not appear to change significantly 
up to 1.5T C [221. l23l. l24l. [25l. [26| in agreement with the small variation of M(T) in the potential model using the Wi(r, T) 
potential as shown in Table I. The widths of many color-singlet heavy quarkonia are broadened suddenly at various 
temperatures [22|, [23|, [24|, HH, [26j|. From the shape of the spectral functions, the range of temperatures from 1.62T C 
to 1.70T C , in which the J/ip width is broadened suddenly, corresponds to the range of J/tjj spontaneous dissociation 
temperatures. Dissociation temperatures for \c and \b in spectral analyses in quenched QCD have also been obtained 
[13, HE HH ■ Dissociation temperatures obtained with the W±(r,T) potential of Eq. (|58|) as well as the F\(t,T) and 
Ui(t,T) potentials are given in Table II. The dissociation temperatures obtained in the Wi(r,T) potential are found 
to give the best agreement with the dissociation temperatures obtained in lattice gauge spectral function analyses, 
as shown in Table II [ll|, [l2| . It is therefore of great interest to see whether the correlators obtained from such a 
potential agree with those from lattice gauge calculations. 

VII. EVALUATION OF THE MESON CORRELATOR IN THE POTENTIAL MODEL 

Upon limiting our attention to the pseudoscalar charmonium with L=0, we note that the quark-antiquark potential 
itself does not possess a potential barrier, and thus there arc no S'-wave resonances [4l[. We calculate the S'-wave 
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phase shift using the amplitude-phase method of Wheeler [47| and Calogero [45|, |48j. The 5-wave phase shifts as a 
function of the continuum energy e and temperature T/T c are shown in Fig. 2. We note that the phase shifts behave 
in two different ways depending on whether there arc bound S states. 




0.5 1.0 

e (GeV) 



FIG. 2: Phase shift as a function of the composite continuum energy e and the plasma temperature calculated in the potential 
model using the Wi(r, T) potential. 



When there are m bound S states in the potential, the S wave phase shift at zero energy will start at rrnr, in 
accordance with the Levinson theorem. For the W\(r,T) potential, there is one bound S-wave state for temperatures 
below ~1.6T C , as indicated by the phase shift of S(e = 0)=ir. The phase shift gradually decreases as energy increases 
and dS(e)/de is negative for all continuum energies. In this case with one or many 5*- wave bound states, the time 
delays for all continuum states are negative and there are no Gamow states in the continuum. Thus, in Eq. (|49|) . only 
a single term, a bound state term with a bound state mass Mb, contributes to the correlator G(r) for T <^ 1.6T C . 
The correlator G(r) of Eq. (|49]) becomes 



G(t) oc 



cosh{M b (r)(r-/3/2)} 
sinh{Af b (T)/?/2)} 



(67) 



where Mb(T) is the mass of the bound S state at the temperature T. A lattice gauge correlator G(r) is represented in 
terms of its ratio with respect to a "reconstructed" correlator G recon (r), which is defined as the correlator calculated 
with the "reconstructed" spectrum at T = 0. In actual practice in the evaluation of the lattice gauge correlators 
GrcconM to obtain the ratios of G(r) /G Tecon (T) shown in Figs, la and 3a, Ref. [24| has used the "reconstructed" 
spectrum at T = 0.75T C which contains only a single bound state [3l|. Therefore, in order to compare with lattice gauge 
G(T)/G Tecon (T), we need to include only a single lowest mass bound state in Eq. (j49|) to evaluate the "reconstructed" 
correlator G Tecon (T) in the potential model. The potential model correlator G(r) relative to the potential model 
G re con(y), normalized to unity at r = 0, is thus 



G(r) 



cosh{M b (T)(r - (3/2)} / cosh{ M b (T = 0)(r - (3/2)} 



Gr, 



i(r) 



sinh{il/ b (T)/3/2)} 



sinh{M 6 (T = 0)(3/2)} 



(68) 



where a = tanh(M b (T)/3/2)/ tanh(M h (T = 0)0/2) - 1. 



Using the above equation that is based on the concept of the absence of L=0 resonance states and Gamow states 
when a bound state is present, we calculate the potential model correlators G(T)/G rccon (T) for the cases of 1.18T C , 
1.25T C , and 1.40T C with bound state masses M{T) obtained in the Wi(r, T) potential as given in Table I. The results 
for these three temperatures are shown in Fig. 3b. As the bound state mass values in Table I have uncertainties 
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FIG. 3: Fig. 3a shows the pseudoscalar charmonium correlators G(r)/G re con(T) obtained in lattice gauge calculations [24 ]. 
The potential model pseudoscalar charmonium correlators G(r)/G re con(T) obtained with the Wi(r,T) potential of Eq. (|58|l are 
shown in Fig. 3b, with the free energy f\(r, T) potential in Fig. 3c, and with the internal energy energy Ui(r, T) potential in 
Fig. 3d. 



of about ±0.7% in M(T), the uncertainties of the potential model correlators arc about ±0.007 in G/G lccon . The 
potential model correlators can be compared with the correlators obtained in lattice gauge calculations show in Fig. 



3a [24[. The general features of the potential model correlators at these three temperatures below about 1.6T C agree 



with those of the lattice gauge correlators. In particular, the slopes d[G /G vccon \/ dr in Fig. 3b for T <~ 1.6T C in 
potential model calculations are of the order ~ 0.03 — 0.1/fm, which do not differ much from the nearly zero slopes 
of d[G/G recon ]/dT in Fig. 3a for T = 0.9T C and 1.5T C in lattice gauge calculations. They differ significantly from the 
general shapes of the potential model correlators obtained by the authors of Ref. 0, [2l[ in Figs, lb and lc, where 
d[G/G rccon ]/dr — 1/fm in Fig. lb and ~ 2.5/fm in Fig. lc, for t < 0.2fm. 

As the temperature increases above ~1.6T C , the 5-wave state calculated with the Wi(r, T) potential is no longer 
bound and the phase shifts at T/T c = 1.95 and 2.60 are shown in Fig. 2. For these cases without a bound 5-wave 
state, the phase shift starts at zero at zero continuum energy e and it increases rapidly as the energy increases. After 
reaching a peak value below 7r/2, the phase shift decreases slowly as the energy increases. Thus, there is a region of 
continuum states for which 38(e) /de is positive. They possess positive time delays and lifetimes. They are Gamow 
states capable of propagating as a composite meson to contribute to the meson correlator. In this case without an 
S-wave bound state, there are no contributions from bound states to the correlator G(t). Only continuum Gamow 
states represented by the integral over k contribute to the meson correlator G(r) in Eq. (jl9")l . 
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0.000 0.010 0.020 0.030 0.040 

e (GeV) 

FIG. 4: (Fig. 4a) The phase shift 8(e) as a function of the continuum energy e at T = 1.95T C and 2.6T C , calculated in the 
Wi(v,T) potential of Eq. ([58}. (Fig. 4b) The width of the Gamow states T(e) as a function of the continuum energy e. (Fig. 
4c) K-(e) as a function of e. 

An expanded view of the phase shift 8(e) as a function of the continuum state k with energy e = k 2 /2/i m is shown 
in Fig. 4a for T = 1.95T C and T = 2.60T C . The corresponding width F(e) extracted from the time delay 88(e) /de) is 
shown in 4b. At T = 1.95T C , the width is zero at e=0, and the gradual increase in width turns into a rapid increase as 
e increases. The increase in width is less rapid at the higher temperature of 2.6T C . Above the energy e(max) at which 
r(e(max)) = r max where r max = 2T, the width of the continuum state is either too large or 88(e) /de is negative, and 
continuum states with energy above e(max) will not contribute to the meson correlator in the integral over k. For 
T = 1.95T C , e(max) is 0.018GeV, and for T = 2.60T C , e(max) is 0.033GeV. They are located at energies slightly above 
the continuum threshold. 

We can also evaluate the amplitude of the continuum wave function at the origin. The comparison of the absolute 
square of the amplitude at the origin r = 0, with the absolute square of the amplitude at very large r gives the K- 
factor for the continuum state [45| shown in Fig. 4c. This /C-factor is quite large at low energies, due to the attractive 
interaction between the quark and the antiquark. The /C-factor decreases as a function of the continuum energy e. 
We need to specify the width limit r max (T) = l/r min lifc (T) = 2T by using Eq. (gSJ). For T = 1.95T C , r min ufe is 
0.185fm and r max is 1.05GcV. For T = 2.6T C , T rain lifc is 0.141fm and T max is 1.40GeV. 

After setting the limit on the width of the Gamow states, we can calculate the pseudoscalar meson correlator, 
G(t) I G vecon (T) normalized to unity at r = 0, for the case without a bound state. The results of the potential model 
correlators for 1.95T C and 2.6T C using the W\(r, T) potential are shown in Fig. 3b. As one observes, the general trends 
of the correlators at these temperatures agree with those from the lattice gauge correlators. 

The results in Figs. 3a and 3b indicate that when the contributions from bound states and continuum states are 
properly treated, there is indeed agreement between the lattice gauge correlators and the potential model correlators 
using the W\(r,T) potential that is a linear combination of F\(r,T) and U\(r, T). 

It is of great interest to investigate whether the potential model correlators depends on thepotential. Accordingly, 
we evaluate the potential model correlators using the color-singlet free energy f \ (r, T) @, S, E3 an d the internal 
energy U\(r,T) 0, [l(| [H, HH for comparison. We need the bound state masses at various temperatures to evaluate 
the correlators, as required by Eqs. (|49|) and (|68|) . We list state energy e(T) and bound state mass M(T) calculated 
in the F\(r,T) potential and the U\(r,T) potential for the L—0 charmonium in Table III. 

Table III. Bound state energy e(T) and bound state mass M(T) (all in GeV) for the L=0 charmonium 
calculated with the F\ and U\ potentials in quenched QCD. 
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T 


1.13T C 


1.18T C 


1.25T C 


1.4T C 


1.6T C 


1.95T C 


2.60T C 


F\ Potential 


e(T) 


-0.0117 


-0.0051 


-0.0018 


-0.000007 


unbound 


unbound 


unbound 


Mass M (T) 


3.029 


3.012 


2.997 


2.972 


U\ Potential 


e(T) 


-0.7483 


-0.4422 


-0.2613 


-0.1331 


-0.07820 


-0.0669 


-0.0225 


Mass M (T) 


3.188 


3.279 


3.285 


3.279 


3.265 


3.259 


3.201 



We note that as a function of temperature T, the L = charmonium masses with the F\ potential are lower than 
the mass at T = and they decrease slowly with temperature up to its dissociation temperature at ~ 1.4T C . On the 
other hand, for the U\ potential the masses are greater than the charmonium mass at T = and vary only slightly 
with temperature. 

We need to calculate further the phase shifts and the wave function amplitude at the origin and evaluate the 
correlators in these potential models. The potential model correlators obtained with F\ and U\ as the potentials are 
shown in Fig. 3c and 3d respectively. The comparison of these potential model correlators with the lattice gauge 
correlators indicates that the correlators calculated with the F\{r,T) and Ci(r,T) potentials deviate from the lattice 
gauge correlators, the deviation being greater for the U\ potential than the F\ potential. The comparison shows that 
among the three different potentials, correlators obtained with the W\(r,T) potential of [O, E2, 0, [3 give the best 
agreement with the lattice gauge correlators. 

VIII. IMPLICATIONS ON THE LATTICE GAUGE SPECTRAL FUNCTION ANALYSIS 

The analysis of spectral function at finite temperatures in lattice gauge theory contains systematic uncertainties and 
lattice artifacts. In the maximum entropy method used for the spectral function analysis, it is necessary to assume 
a default spectrum to define the entropy. The spectral function obtained in the maximum entropy method depends 
on the assumed default spectrum (see Fig. 4 of [2J|)- The extracted spectra often exhibit two prominent broad peaks 
in the continuum which do not seem to correspond to physical continuum states 0, HH, 0, [25|, [26[. While it is 
important to study theoretically what the lattice artifacts are expected to be, some prior knowledge of the spectral 
function as inferred from the potential model will be useful to provide additional information for the assumed default 
spectrum and the expected spectral function. 

From the physical picture that emerges in the evaluation of the current-current correlator, the spectral function 
analysis in the lattice gauge theory should be guided by the basic principle that out of the set of continuum states, only 
resonance states and Gamow states with lifetimes of sufficient magnitudes can propagate as composite objects and 
can contribute to the current-current correlator. The potential model can provide useful information on the nature 
of bound states, resonance states, and Gamow states for spectral function analysis. We can discuss first the case of 
L = 0. As the Q-Q interaction itself does not possess a barrier, we expect that there are no L = potential resonances 
in the continuum [4l| . The occurrence of L = bound states in the potential is accompanied by phase shifts with 
negative 85(e) /de (Fig. 2) with no delay time and composite lifetime for the continuum states. Consequently when 
L = bound states occur, only bound states can contribute to the correlator. According to these results of the 
potential model, the lattice gauge spectral function analysis for L = states should include only bound states without 
continuum states. On the other hand, when the L = bound states are absent, states with delay time and widths 
occurs at energies close to e ~ 0. The correlator receives contribution only from this region of continuum states (see 
Fig. 3b) and lattice gauge spectrum function analysis for L = states should employ a default spectrum that contains 
this element of the continuum contribution. 

For the case of L > 0, the presence of centrifugal barrier leads to potential pockets and possible potential resonance 
states, which may or may not occur, depending on the strength of the potential and the centrifugal barrier. One 
expects that when potential resonance states occur, the phase shift decreases as a function of e at continuum energies 
above the resonance energies, and only bound states and possible resonance states contribute to the correlator in 
lattice gauge analysis. On the other hand, when bound states and resonance states are absent (as for example at 
high temperatures), Gamow states with delay time and widths of various magnitude occurs. The correlator receives 
contribution only from these Gamow states. Lattice gauge spectrum function analysis for L > states should employ 
a default spectrum that contains these essential characteristics. 

IX. CONCLUSIONS AND DISCUSSIONS 

The potential model has been a useful concept in the physics of heavy quarkonia since the discovery of J/ip. It 
provides a useful tool to examine quarkonium energies, quarkonium wave functions, reaction rates, transition rates, 
and decay widths [l], 0, H, 0, H[ . It allows the extrapolation to the region of high temperatures by expressing screening 



16 



effects in terms of the temperature dependence of the potential @, 0, S, U, EE EL EE EE El EE EE EE EE El ■ The 

comparisons of the potential model dissociation temperatures with the lattice gauge spectral function dissociation 
temperatures show consistency when one uses the potential that is a linear combination of the free energy and the 



In addition to lattice gauge spectral function analyses, results have also been obtained for the lattice gauge current- 
current correlators [20I . l2lj . The current-current correlator is related to the spectral function by a generalized Laplace 
transform. In principle, they carry equivalent information. One expects that the consistency of the potential model 
with lattice gauge spectral function analyses should be extended to the comparison of thepotential model current- 
current correlators with lattice gauge current-current correlators. Recent works in Refs. [20I . l2lj however make the 
contrary claim that the meson correlators obtained from many different types of potential models arc not consistent 
with lattice gauge correlators and potential models cannot describe heavy quarkonia above T c . In the work of Ref. 
pol . [2l| . continuum states arising from a free fermion Q and Q pair in a fermion gas contribute to the correlator. 
However, to be able to propagate as a composite meson so as to contribute to the correlator, the quark and antiquark 
must be temporally and spatially correlated. Continuum states in the free quark and antiquark gas may not have 
sufficient temporal and spatial correlations to qualify as composite objects for such a propagation. 

Following the basic principle that among the continuum states only resonance states and Gamow states with a 
lifetime of sufficient magnitude can propagate as a composite object and contribute to the correlator, we re-evaluate 
the current-current correlator in the potential model. With the simple example of the pseudoscalar correlator, we 
show in this paper that when the bound state and continuum state contributions are properly treated, the W\(r,T) 
potential using a linear combination of F\ and U\ proposed in [TTI.[T3| gives correlators consistent with those of lattice 
gauge correlators, while the F\ potential and the U\ potential lead to deviations. Our results indicate consistency 
of the W\(t,T) potential proposed in [HI El with both the lattice gauge spectral function analysis and the lattice 
gauge correlator analysis. The present agreement is not surprising as the current-current correlator and the spectral 
function are related by a generalized Laplace transform, and they indeed carry equivalent information. 

There are uncertainties, limitations, and future prospects in the potential model that requires further investigations. 
Lattice gauge calculations provide information on thermodynamical quantities of the free energy F\ and the internal 
energy U\. We have shown deductively in Eq. (11) of Ref. [llj that the internal energy U x {r,T) contains U^(r,T) 

and Ug(r,T) — U g a(T). Only Uqq(t,T) pertains to the interaction potential between Q and Q. The other parts 
need to be subtracted out from U\(r,T) to obtain the Q-Q potential. We have made use of the knowledge of the 
equation of state from an independent lattice gauge calculations to evaluate U^l.(r,T) from C/i(r,T), leading to the 
present linear-combination model of the U^L(r,T) potential of Eq. Q58|) . However, a more rigorous treatment will 

involve a direct evaluation of the Q-Q potential by evaluating the quantities of Ug (r, T) — U g o(T) in the lattice gauge 
theory. It will be of great interest to examine how one can determine directly the Q-Q potential in a lattice gauge 
calculations without resort to the equation of state. 

The potential model we have developed so far has the limitation that the important spin-spin interaction has not 
been included. As the spin-spin interaction is responsible for the r] c -J '/ip and rj^-T splittings, it modifies the binding 
energies and the dissociation temperatures. It is important to include spin-spin and other components of the Q-Q 
interactions in the potential model to see how they may affect the stability of heavy quarkonia. 

The potential model allows the evaluation of many quantities of interest. The potential model in the present 
manuscript uses thermodynamical quantities obtained in quenched QCD. The effects of dynamical quarks on the 
stability of quarkonia can be studied by using potential models extracted from thermodynamic quantities calculated 
in full QCD EE HH. One can examine the quark mass dependence of quarkonium stability and quarkonium gluon 
dissociation cross sections, some results of which have been presented recently EE E3|- The heavy quark potential 
model so far developed has the limitation that it is restricted to non-relativistic cases. It will also be of interest to 
study the relativistic effects by examining relativistic two-body potential models along the lines of Dirac's constraint 
dynamics as in Ref. |,1 , which will allow us to study the stability of light-quark systems within the potential model. 

After the present manuscript has been submitted for publication, a recent preprint [il refers to our procedure of 
including only states of sufficient lifetimes and asserts that "this procedure might be incorrect, since the evaluated 
correlator has to be compared with the lattice ones, which do have a free gas (infinite temperature) limit." 

This statement in Ref. [il with regard to our work might be incorrect. Firstly, the evaluated lattice gauge 
correlators of Ref. [11 to be compared were carried out in finite temperatures with T ~ (1 — 3)T C and T « r M 
(Too —* 00) • This temperature T is in the non-perturbative region and is far from the infinite temperature perturbative 
QCD limit. The spectrum of a free gas pQCD limit at infinite temperature is not relevant to the finite temperature 
lattice gauge correlators being compared at hand. Secondly, the presence of a continuum spectrum in the infinite 
temperature limit is in fact consistent with our procedure of including states of sufficient lifetimes at that infinite 
temperature limit. 
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The arguments to support our procedure have been presented already in the manuscript. We repeat the main points 
again below to rebut the statement of (49j . 

As discussed in Eq. (1), the meson correlator describes the probability amplitude for creating a composite meson 
at time and subsequently destroying the composite meson at time r. The operation of destroying the meson at time 
t can be considered as an operation of measurement (or an operation of detection) of the meson at time r. From 
the discussions in Eqs. (28) and (29), r is less than and up to (3/2 or 1/2T, where T is the temperature, and the 
measurement operation takes place within r < 1/2T. At the temperature T, in order to be detected by the correlator 
measurement at r = 1/2T, the composite meson state needs to have a lifetime Tiif e greater than 1/2T, which has been 
conveniently labeled as r min iif e (T) = 1/2T in Eq. (|45|) . the minimum meson state lifetime (for correlator detection) 
at temperature T. 

From this analysis, the composite meson states that can be detected by the correlator measurement will include 
meson states of shorter and shorter lifetimes as temperature increases, and the correlator spectrum above the bound 
states region will depend on the temperature. In the infinite temperature limit T = Too — > oo, the minimum lifetime 
7"min life (Too) ? which is equal to 1/2TQO, approaches zero. The correlator measurement includes states that live for a 
very short composite object lifetime th[ c > T m - U1 iif c (Too) ~ 0. This condition can be satisfied for weakly interacting 
free gas continuum states. The correlator allows a free gas continuum spectrum in the infinite temperature limit. 
The presence of a continuum spectrum in the infinite temperature limit is consistent with our procedure of including 
states of sufficient lifetimes. 

We turn our attention now to the finite-temperature lattice gauge calculations in Ref. [24] with T ~ (1 — 3)T C 
and T « Too- Our afore-mentioned comparison of the meson lifetime and the minimum lifetime for the correlator 
measurement indicates that only composite states that have lifetimes rnf e greater than r m j n nf c (T) = 1/2T can survive 
and be detected by the correlator measurement at T. Free gas continuum states with very short composite object 
lifetimes T m - ln iif e (Too) = 1/2X00, much smaller than T m ; n iif c (T) = 1/2T, cannot survive up to r m ; n m c (T) = 1/2T and 
will not be detected by the correlator measurement at T. Because of this limitation, the meson spectrum obtained 
in finite-temperature correlator calculations at T should be different from the correlators in the infinite temperature 
limit and should only contain bound states, resonance states, and Gamow states of sufficient lifetimes greater than 
1/2T. 

Based on the above rebuttal, there might be no basis for the authors in Ref. [49| to make the statement mentioned 
above. 

In conclusion, we have found that the potential model of Ref. I 111 is consistent with both spectral function analyses 
and current-current correlator analyses. The potential model of [111 ] can be a useful tool to complement lattice gauge 
calculations for the study of heavy quarkonia at high temperatures. 
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